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The orbital motion of the Laser Interferometer Space Antenna (LISA) introduces modulations 
into the observed gravitational wave signal. These modulations can be used to determine the 
location and orientation of a gravitational wave source. The complete LISA response to an arbitrary 
gravitational wave is derived using a coordinate free approach in the transverse-traceless gauge. The 
general response function reduces to that found by Cutler U for low frequency, monochromatic plane 
waves. Estimates of the noise in the detector are found to be complicated by the time variation of 
the interferometer arm lengths. 



I. INTRODUCTION 

In this paper we derive the response of the Laser Inter- 
ferometer Space Antenna (LISA) to an arbitrary gravi- 
tational wave. Our detector response function includes 
the full orbital motion and is valid at all frequencies. 
Previous treatments have either assumed that LISA is 
stationary with respect to the background sky0j0, Q , or 
have been limited to the low frequency limit [ij, |5j • 

The LISA mission^ calls for three spacecraft to or- 
bit the Sun in an equilateral triangular formation. The 
center-of-mass for the constellation, known as the guid- 
ing center, follows a circular orbit at 1 AU and has an 
orbital period of one year. In addition to the bulk motion 
of the detector about the Sun, the triangular formation 
will cartwheel about the guiding center in a clockwise 
manner as seen by an observer at the Sun. The orbital 
motion introduces frequency (Doppler), amplitude, and 
phase modulation into the observed gravitational wave 
signal. These effects have been calculated in the low fre- 
quency limit, where the antenna pattern is well approx- 
imated by a quadrupole, and the Doppler modulation is 
due to the guiding center motion. At higher frequencies 
the antenna pattern becomes more complicated, and the 
"rolling" Doppler modulation due to the cartwheel mo- 
tion has to be included. 

The divide between high and low frequencies roughly 
coincides with the divide between gravitational waves 
with wavelengths shorter or longer than the arms of 
the detector. To be precise, the dividing line is the 
transfer frequency /, = c/{2ttL), where L is the unper- 
turbed distance between spacecraft (3| . For LISA, with 
its mean arm length of 5 x 10^ km, the transfer fre- 
quency has an approximate value of 10 mHz. Above 
the transfer frequency the antenna pattern is distinctly 
non-quadrapolar0. The frequencies at which the guid- 
ing center and rolling motion impart measurable effects 
are easily estimated. Since both the guiding center and 
cartwheel motions have periods of one year, the Doppler 
modulations will enter as sidebands separated by the 
modulation frequency /,„ = 1/year. Equating the mod- 
ulation frequency to the Doppler shift, Sf ~ {v/c)f, for 
motion with velocity v yields the characteristic frequency 
fv — cfm/v at which Doppler modulation becomes mea- 
surable. Assuming 5 x 10^ km armlengths the cartwheel 



turns with velocity v/c — 0.192 x 10^^, while the guiding 
center moves with velocity v/c = 0.994 x lO^**. Thus, 
the Doppler modulation due to the guiding center's mo- 
tion becomes measurable at frequencies above fgc = 0.3 
mHz, while the rolling cartwheel motion becomes impor- 
tant above fr = 16 mHz. This demonstrates that at low 
frequencies only the bulk motion of the detector needs 
to be considered, while at high frequencies the cartwheel 
motion also needs to be included. 

Our calculations are performed using the coordinate 
free approach introduced in Ref.0| (see also Ref. for 
a closely related approach) . This allows our results to be 
applied to any variation on the current LISA design, or 
to any follow on mission. The low frequency limit of our 
general detector response function yields a simple result 
that can be shown to agree with Cutler'sT|. Throughout 
this paper we use natural units where G — c — 1 , however 
we will report all frequencies in terms of Hertz. 



II. DETECTOR RESPONSE 

The detector response to a gravitational wave source 
located in the n direction can be found using Barycentric 
coordinates {t, x) and the transverse-traceless gauge to 
describe a plane gravitational wave h(g, il) propagating 
in the 51 = — n direction. The surfaces of constant phase 
are given by ^ = t+fi-'x. = const.. A general gravitational 
wave can be decomposed into two polarization states: 

h(e,n) =/^+(Oe+ + /ix(e)e" (1) 
where e+ and arc the polarization tensors 

e+ = ii ^ u — V ^ i), 

= M (g) w + u (g) M. (2) 

The basis vectors u, v and the source location n can be 
expressed in terms of the location of the source {6, (p) 
according to 

u — cos 6 cos (f>x + cos 9 sin (t>y — sin 9 z 
V = sin 4>x — cos (p y 

n = sin 9 cos 4>x + sin 9 sin (py + cos 9 z . (3) 
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Following the Doppler tracking calculations described 
in Ref.Q, we find that the optical path length between 
spacecraft i and spacecraft j may be written as 



= \^jitj)-Mu)\ 



+ ^ihAU)^hjiU)): J h(e(A))dA. (4) 

Here = 1]^,^ + h^^ and |x — y| denotes the Cartesian 
distance between x and y. The unit vector 



(5) 



points from spacecraft i at the time of emission, U, to 
spacecraft j at the time of reception, tj. Finally, the 
quantity ^(A) is the parameterized wave variable 



C(A) = i(A) -f7-x(A), 



(6) 



Explicitly, the time and position depend on the parame- 
terization in the following way 



t(A) ^ U + X 
x(A) = Xi{U) + \fij{ti) . 



(7) 

(8) 



The variation in the Cartesian distance between the 
spacecraft can be separated into a contribution due to 
the Sun and other Solar system bodies, and a small per- 
turbation due to the gravitational wave. Denoting the 
unperturbed spacecraft locations by x°(i) and integrat- 
ing the geodesic equation for the metric g^j/ yields 



|x,(t,)-x,(t,)|-K(t,)-x°(t,)| 



dt 



dxj (t) 

-h, iu) ■ J\{t~n- x,(t)) • dt . (9) 

The gravitational wave dependent terms in the above 
equation can be ignored as they are down by a factor 
of u = \dx{t)/dt\ ~ 10^'' compared to the leading order 
gravitational wave contribution described in Q). Thus 
we may write 



where 



£°,(t.) = |x°fe)-x^(t.)|. 



and 



1 fij{ti) (g) fjjiU) 

2 i-n-nj{t,) 



HOd^. 



(10) 



(11) 



(12) 



Here we have used 

dX 



= 1 - • fijiu) 



(13) 



to make a change of variables in the integration. Defining 
H(a,&) to be the antiderivative of the gravitational wave 



H(a,6) 



simplifies our expression further to 

1 r,j{t,) (g) f,j{t,) 



2 l + n-f,j{ti) 



(14) 



(15) 



where we have used the relationship Q — —n. To leading 
order in h, the time of reception tj is defined in terms of 
the time of emission ti by the implicit relation 



£0(i,) = |x°(t.+e(iO)-x°(tOI 



(16) 



The gravitational wave can be decomposed into fre- 
quency components 



h(^) = / h{u;)e"^^duj 

J —OO 

which allows us to write 

h(c.) 



H(a,6) 



(17) 



) du (18) 



and 



5e^J{ti) = £^j{t^) / D(w,t„n) : h{iu)e''^^' dcu (19) 



(20) 



where the one-arm detector tensor is given by 

T>{u;,ti,n) = i {fij{t,) (E)f,j{ti))T{uj,t„n) 
and the transfer function takes the form 

UJ 



T{uj,t^,h) 



smc 



X exp 



2uj. 



-(1 -I- n • fij{ti)) 



2wii 



(1 -I- n • hjiti)) 



(21) 



Here = l/iij{ti) is the angular transfer frequency for 
the arm. 

The connection between the optical path length vari- 
ations and the detector output depends on the interfer- 
ometer design. The original proposal was to use laser 
transponders at the end-stations to send back a phased 
locked signal. A more recent proposal is to eliminate the 
transponders and turn LISA into a virtual interferometer 
where the signal is put together in software. The raw in- 
gredients for this procedure are the phase differences be- 
tween the received and transmitted laser light along each 
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arm. The signal transmitted from spacecraft i that is re- 
ceived at spacecraft j at time tj has its phase compared to 
the local reference to give the output ^ij{tj). The phase 
difference has contributions from the laser phase noise, 
C{t), optical path length variations, shot noise and 
acceleration noise n" (t) : 

+ntj [t, ) - iU) ■ {K^ {t, ) - n% iU)) . (22) 

Here U is given implicitly by U — tj — iij{ti) and vq is 
the laser frequency. We have included the variations in 
the optical path length caused by gravitational waves, 
SiijiU), and those caused by orbital effects, A£ij{ti). In 
what follows we will ignore the orbital contributions to 
the phase shift as they can be removed by high pass fil- 
tering. The subscripts on the noise sources identify the 
particular component that is responsible: Ci is the phase 
noise introduced by the laser on spacecraft i, nfj denotes 
the shot noise in the photodetector on spacecraft j used 
to measure the phase of the signal from spacecraft i, and 
denotes the noise introduced by the accelerometers 
on spacecraft j that are mounted on the optical assem- 
bly that points toward spacecraft i. 

The three LISA spacecraft will report six phase differ- 
ence measurements which can then be used to construct 
a variety of interferometer outputs. The simplest are the 
three Michelson signals that can be formed by choosing 
one of the spacecraft as the vertex and using the other 
two as end-stations. The Michelson signal extracted from 
vertex 1 has the form 

Slit) = $12(^2) + *2i(t) - $13(^3) - $3i(t) , (23) 
where t2 and ^3 are given implicitly by 



t2 



(24) 



Unfortunately, the Michelson signals will be swamped by 
laser phase noise, so a more complicated virtual interfer- 
ometer signal has to be used. The X variables are a set 
of three Michelson-like signals that cancel the laser phase 
noise .^j. The X signal extracted from vertex 1 has the 
form 

Xi{t) ^ $12(t2)- $12(^1+^2-0 
+ $2l(i) - *2l(tl) 

- $13(^3) + *13(il+t3-0 

- $3i(i) + $3i(ii). (25) 



Given a gravitational wave signal h(g, fi), a model of the 
instrument noise, and a description of the interferome- 
ter's orbit, we can use equations (|22|l and (|25|l to 
calculate the detector response. Since the entire calcu- 
lation is performed in Barycentric coordinates, the time 
t that appears in H23|) and H25|) is not the time r mea- 
sured by the clock on spacecraft 1. They are related by 
the standard time dilation formula dr = dt^Jl — v\{t). 
However, since we only need to work to leading order in 
w, there is no need to distinguish between t and r. 

Our expression for the LISA response is much more 
complicated than the previous approximate descriptions. 
The time- variation of the optical path lengths is the main 
cause of the difficulty. It is responsible for the implicit 
relations that riddle the calculation. The path length 
variations have three main causes - intrinsic, tidal, and 
pointing. The intrinsic variations are part and parcel of 
the cartwheel orbit, which only keeps the distance be- 
tween the spacecraft constant to leading order in the or- 
bital eccentricity. The tidal variations are caused by the 
gravitational pull of other solar system bodies, most no- 
tably the Earth and Jupiter. The pointing corrections 
are a relativistic effect caused by the finite propagation 
speed of the lasers, which means that the spacecraft move 
between transmission and reception. The latter effect can 
be separated from the others: 

^^3{t^) = L,j{ti) [1 + hj{ti) ■ V,(t,) + Oiv"^)] (27) 

where Vj is the velocity of spacecraft j and 

Uj{t,) = |x,(t,)-x,(t,)|. (28) 

Ignoring tidal distortions and working to second order 
in the orbital eccentricity e, the orbits described in the 
Appendix yield 



Li2{t) = L [I + 15sin(a-f ^) -cos(3a) 



(29) 



and similar, yet slightly different, expressions for Li^{t) 
and £23 (i)- Here a{t) — 2tt f„it 4- k is the orbital phase 
and fm — 1/year is the orbital frequency. The mean 
armlength L is related to the eccentricity e and semi- 
major axis a by L = 2-\/3ae. Setting L = 5 x 10^ meters 
yields e = 0.00965 « 10~^. The spacecraft have velocities 
of order v « 2nfma ~ 10^^. Using these numbers we see 
that the lowest order intrinsic variation is far larger than 
the pointing variation. The tidal variations turn out to 
be comparable to the intrinsic variation 6J and therefore 
should not be ignored. 



III. STATIC LIMIT 



The times ti, t2 and are defined implicitly: 

h = t - ei2{ti) - e2i{t2) 

= t - luiti) - isiits) 

t2 = t~ £21(12) 

ts = t-£3l{t3). 



(26) 



As a point of reference, we can apply our general 
method to a static, equal arm detector interacting with a 
monochromatic, plane-fronted gravitational wave propa- 
gating in the f2 — —h direction with principle polariza- 
tion axes p and q: 



h(/,0- A 



A 



g27ri/(t-|-n-x)g> 



(30) 
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where and are complex constants and 



p® p ~ q® q, 
p® q + q®p. 



(31) 



Defining tlie polarization angle ■0 — — arctan('C • p/u • p) 
we have 



cos 2^1) e"*" — sin 2i/) 
sin 2-0 + cos 20 



(32) 



Thus, in terms of the general decomposition Q we have 

h+{i) = (A+cos20 + Axsin20)e2"/(*+"-^), 
hAt) = (AxCOs20-^+sin20)e2'^*^(*+"-''). (33) 

The signal portion of the X variable defined in H25|) 
reduces to 

Xl{t) = 2Tiv^{5ti2{t - 2L) - Sii2it - 4L) 
+Si2i{t - L) - Sl2i{t ~ 3L) ~ S£i3{t - 2L) 
+Sh3it - 4L) - Si3i{t -L)+ Stii (t - 3L)) (34) 

In terms of the strain xi{t) ~ Xf{t)/{2TTLh'o) we have 



t:i(t) = D(n,/) :h(/,0sin2(///,) 



(35) 



where 



D(n,/) 



1 



((^12 «> ?'l2)^(^12 • n,f) 
(ri3 «) ^13) r(f 13 -n, /)) 



(36) 



and 



m/) 



/(1 + s) 
smc I — — ) exp 



-I — (3 



, 7(l-^) \ ( ■ f n , 
smc cxp —I 1 — s) 

' 2/, ; 2// ' 



(37) 



Orienting the detector in the x — y plane according to 
Figure 2 of Ref . Q , we have 

fi2 = cos(7r/12) X + sin(7r/12) y 

fi3 = cos(57r/12) X + sin(57r/12) y . (38) 

Combing these expression with equations Q and Q 

yields 



(''12 ®ri2) : 
{fi2 ® fi2) : e' 
(^13 ® ris) ■ 



1 



((1 + cos2 9) sin(20 + 7r/3) — sin 0] 



cos6 sin(20 — 7r/6) 

- ((l + cos^ 6) sin(20-7r/3) — sin 9) 



(fi3 (g) fi3) : = — cos6' sin(20 + 7r/6). 



(39) 



The above collection of equations, (^5)1 through (pT^ . 
fully define the detector response in the static limit. The 



X variable response xi (t) is related to the Michelson re- 
sponse Si{t) by 



(t) = x,{t) sin-2(///,) = D(n, /) : h(/,0 . (40) 



Our expression H4U|) for si{t) agrees with that quoted in 
Ref. 4]. In the low frequency limit, f f*, the transfer 
function reduces to unity, T — I, and 

si(t) = (A+F+(0,0,^)+Axi^^(0,0,0'))e2-^*. (41) 

The beam pattern factors 

F+(6I,0,V) = ^(?'i2®ri2-fi3«)fi3) :e+ 
F^(0,0,V') = ^(ri2®ri2-fi3®fi3) (42) 

take their familiar form: ^3 

= ^fi(l + 008^6') cos 20COS2V' 



2 V2 



COS 9 sin 20 sin 2-0 



/3 1 

F"" = ^f-(l + cos2 6')cos20sin20 



2 V2 



cos 9 sin 20 cos 2iIj 



(43) 



Notice that the overall factor of ^/S/2 compared to a 
detector with 90° arms come out naturally in our calcu- 
lation. 



IV. LOW FREQUENCY LIMIT 

As a further point of reference, we can apply our 
general result to the low frequency limit considered in 
Ref. l] . The low frequency limit is defined by the condi- 
tion / ^ /*, where /* = l/(27rL) is the typical trans- 
fer frequency along each arm. A LISA mission with 
L = 5 X 10^ meter arms has a transfer frequency of 
/, = 0.00954 « 10-2 Hz. The motion of the LISA con- 
stellation is included to leading order in the eccentricity, 
and the gravitational wave is taken to be monochromatic, 
plane-fronted and propagating in the = —n direction: 

h(t,x) = A+e+cos(27r/(t+ri-x))+Axe'' sin(27r/(t+n-x)) . 

(44) 

The two orthogonal polarizations have constant real am- 
plitudes A+ and Ax ■ The basic Michelson signal consid- 
ered by CutlerU takes the form 



Slit) = 



Sii2{t - 2L) + S£2i{t - L) 
2L 

Shsjt - 2L) + Sl^ijt - L) 
2L 



(45) 



This expression ignores the time variation of the arm- 
lengths due to higher order terms in the orbital eccen- 
tricity or perturbations from other solar system bodies. 
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Using lfTK|l we find 



and 



si{t) = F+{t)A+coa{2Trf[t + n-xi{t)]) 

+ ^^^(i)AxSin(2^/[t + n-xi(t)]), (46) 



where 



F+{t) = - {cos2^l; D+{t)~ sin D'^it)) 

= ^ {sin 211; D+{t) +008 24! D'^it)) , (47) 



and 



D+{t) = {fi2{t) <E) h2{t) - hsit) <E> hsit)) : e+ 
D'{t) = (fi2(t)®fi2(t)-fi3(t)®fi3(i)) :e\(48) 

The expression for the strain in the detector can be re- 
arranged using double angle identities to read: 

si (t) = A{t) cos [2nft + [t) + dpp {t)] . (49) 

The amplitude modulation A{t), frequency modulation 
ipoit) and phase modulation 0p(i) are given by 



Ait) = [{A+F+it)f + {A,F^t)) 

hoit) - 27r/n-xi(i) 

— 211 fa sin 6 cos(a — 0) 



21 1/2 



(j)p{t) = — arctan 



Axf"(t) 

A+F+{t) 



(50) 



(51) 



(52) 



Using the orbits described in the Appendix, the coordi- 
nates of each spacecraft are given to leading order in the 
eccentricity by 



X = a cos(a) -|- ae (sina cosa sin/3 — (1 + sin^ a) cos/3) 
y — a sin(a) -f ae (sina cos a cos/3 — (1 + cos^ a) sin/3) 
z = — \/3ae cos(a — /3) , (53) 

where a = 2'iTfmt -I- k is the phase of the guiding center 
and /3 = 2n7r/3 -|- A is the relative phase of each space- 
craft in the constellation (n = 0, 1,2). The unit vectors 
fij{t) can be derived from the coordinates given in (|53|l . 
Putting this all together yields 

D+it) = — - 36sin2 6'sin(2a(i) - 2A) 
64 L 

+ (3 + cos 26) cos 20(9 sin 2A - sin(4a(t) - 2A)) 
+ sin 2(/)(cos(4a(t) 2A) - 9 cos 2A)) 

- 4\/3sin26l(sin(3a(t)-2A-0) 

-3sin(a(t) - 2A + 0))j (54) 



I?x(t) 



1 

16 



\/3cos6'(^9cos(2A- 20) 
- cos(4a(t) - 2A - 20)^ 
6sin6'(^cos(3a(i) - 2A - 0) 
+3cos(a(t) - 2A -}-(/))) 



(55) 



Finally, for circular Newtonian binaries, the polarization 
angle tp can be related to the angular momentum orien- 
tation L {Ol,4>l) by 



tan-f/j 



■C • p L ■ u 
u • p L-v 

cos6'cos(0 — 0l) sin^i — cos 9l sin( 



where we have used 



sin 6'l sin(0 - 0l) 



h X L. 



,(56) 



(57) 



The parameters k and A define the initial location and 
orientation of the LISA constellation. They are related 
to the quantities 0o a-nd ag in Cutler'silJ equations (3.3) 
and (3.6) according to k = 0o and A — 37r/4 -I- 0o — 
cko- Our compact expression (|49(l for the low frequency 
limit agrees with Cutler's[3| result, but the agreement 
is by no means obvious. The equality can be established 
using a computer algebra program or by direct numerical 
evaluation. 



V. SPECTRAL NOISE 

The variation in the optical path length will be re- 
flected in the noise transfer functions. For example, the 
noise n|;^ (t) enters into the X variable as 



N{t)^n'^,{t)-n'^,{h). 



(58) 



Writing ?^2i(0 ~ ^(0 and assuming that the armlengths 
are fixed, £31 (ti) = ^31(^3) = L, yields the standard result 



and 



N{f)^n{f) (l-e^W/.^ 



SN{f)^isin^^] Snif). 



(59) 



(60) 



Here S'„(/) is the noise spectral density in the photode- 
tector and we have used the assumption 



<n{f)n*{f')>^Sif-r)SM), 



(61) 



where the brackets <> denote an ensemble average. The 
situation is much more complicated when the armlengths 
vary since 

^(/) = "(/) - ^ / n{h{t))e'-^f'dt . (62) 
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Because £{t) varies with a one year period, the transfer 
function will develop sidebands at / ± n/„ where n takes 
integer values. Working to zeroth order in v and lowest 
order in e we have 

tiit) ~ t - 2Li2{t) . (63) 
Using the expression H29(l for Li2{t) and the expansion 



^ix sin(27T f^nt) 



e ^fcWe'"^^"'*, (64) 



k— — c 



where Jk is a Bessel function of the first kind of order fc, 
allows us to write 



oo oo 



Nif) = nif)- J2 E nif+ij+3k)U)e'^f/f' 

j — — oc k— — oo 

^^2^f^U+3k)/U^^MJ-3k)/6JJ^] ( .f^VfiS) 



16/* 



16/* 



The dependence on the Bessel functions tells us that 
the sidebands only become significant for frequencies ap- 
proaching /*/e ~ 1 Hz. Below the transfer frequency 
/, ~ 10 mHz it is safe to ignore the time variation of the 
arm lengths in calculations of the noise transfer functions. 



and low frequency limits. However, we have said little 
about how the general result should be used. Given a 
specific model for the orbit, such as the simple Keple- 
rian model described in the Appendix, it is possible to 
solve the implicit relations for the detector orientation, 
arm lengths and emission times, as we did in equations 
(|?7jl . and These can then be used to give ex- 

plicit expressions for the Michelson or X variables. We 
did not quote these expressions as they are very large 
and not very informative. Ultimately, any application 
that requires the full LISA response function is likely to 
be numerical. It is a simple matter to write a computer 
program that returns the LISA response function using 
equations (|25|l . H22(l and (|15() . If one is only interested 
in sources with frequencies below ~ 5 mHz, the low fre- 
quency approximation (|49|l will suffice, but for accurate 
astrophysical parameter estimation above 5 niHz, the full 
LISA response function has to be used. 
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Appendix: Keplerian Spacecraft Orbits 

For a constellation of spacecraft in individual Keple- 
rian orbits with an inclination of i = \/3e the coordinates 



of each spacecraft are given by the expressions 
X ~ r ^cos(\/3e) cos/3cos7 — sin/3 sin7 
y = r ^cos(V3e) sin/3cos7 -I- cos/3sin7 
z — — r sin(\/3e) COS7 . 



(66) 



where /3 — 2n7r/3 -I- A (ri = 0, 1, 2) is the relative orbital 
phase of each spacecraft in the constellation, 7 is the 
ecliptic azimuthal angle, and r is the standard Keplerian 
radius 



a(l 



1 + e cos 7 



(67) 



Here a is the semi-major axis of the guiding center and 
has an approximate value of one AU. 

To get the above coordinates as a function of time we 
first note that the azimuthal angle is related to the ec- 
centric anomaly, i/), by 



tan ( ^ ) = 



1 + e /V 

tan — 

1 - e \ 2 



(68) 
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and the eccentric anomaly is related to the orbital phase 
a{t) — 27r/mt + K through 



a ^ P — ip — e sin ^ . 



(69) 



For small eccentricities we can expand equations (|68|) and 
(|69|l in a power series in e to arrive at 



y = asin(Q;) + ae (sinacosacos/3 — (1 + cos^ a) sin/3) 
z = — \/3aecos(a — /3) . (71) 



7 — (a— /3) + 2esin(a — /3) + -e^sin(Q! — /3)cos(a — /?) + . . . 

(70) 

Substituting this series into equation Ht)t)|) and keeping 
terms only up to order e gives us 

X = a cos(a) + ae (sina cosa sin/3 — (1 + sin^ a) cos/3) 



These are the desired coordinates of each spacecraft as 
a function of time. Notice that by keeping only linear 
terms in the eccentricity we are neglecting the variation 
in the optical path length. The path length will change 
due to the Keplerian orbits, but these effects enter at 
0{e^) and above. 



